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Abstract — Non-vascular plants (lichens and mosses) are signifi- 
cant components of tundra landscapes and may respond to climate 
change differently from vascular plants affecting ecosystem carbon 
balance. Remote sensing provides critical tools for monitoring 
plant cover types, as optical signals provide a way to scale from 
plot measurements to regional estimates of biophysical proper- 
ties, for which spatial-temporal patterns may be analyzed. Gas 
exchange measurements were collected for pure patches of key 
vegetation functional types (lichens, mosses, and vascular plants) 
in sedge tundra at Barrow AK. These functional types were found 
to have three significantly different values of light use efficiency 
(LUE) with values of 0.013 ± 0.001, 0.0018 ± 0.0002, and 
0.0012 d= 0.0001 mol C mol -1 absorbed quanta for vascular 
plants, mosses and lichens, respectively. Discriminant analysis 
of the spectra reflectance of these patches identified five spectral 
bands that separated each of these vegetation functional types as 
well as nongreen material (bare soil, standing water, and dead 
leaves). These results were tested along a 100 m transect where 
midsummer spectral reflectance and vegetation coverage were 
measured at one meter intervals. 

Along the transect, area-averaged canopy LUE estimated from 
coverage fractions of the three functional types varied widely, 
even over short distances. The patch-level statistical discriminant 
functions applied to in situ hyperspectral reflectance data collected 
along the transect successfully unmixed cover fractions of the 
vegetation functional types. The unmixing functions, developed 
from the transect data, were applied to 30 m spatial resolution 
Earth Observing-1 Hyperion imaging spectrometer data to ex- 
amine variability in distribution of the vegetation functional types 
for an area near Barrow, AK. Spatial variability of LUE was 
derived from the observed functional type distributions. Across 
this landscape, a fivefold variation in tundra LUE was observed. 
LUE calculated from the functional type cover fractions was 
also correlated to a spectral vegetation index developed to detect 
vegetation chlorophyll content. The concurrence of these alternate 
methods suggest that hyperspectral remote sensing can distinguish 
functionally distinct vegetation types and can be used to develop 
regional estimates of photosynthetic LUE in tundra landscapes. 
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I. Introduction 

H IGH northern latitudes are undergoing dramatic changes 
in climate. Warming trends have been observed in 
northern regions and Global Circulation Model predictions 
indicate arctic and boreal regions are likely to warm by several 
degrees over the next century, a greater rate than other regions 
on the globe [1], [2]. In high latitudes seasonal shifts in surface 
properties involve freezing and thawing of water, so relatively 
small temperature changes around the freezing point can result 
in significant environmental shifts. Along with temperature 
changes, significant changes in precipitation and evapotranspi- 
ration are also predicted for these regions [1], [2]. 

Tundra vegetation is expected to be particularly responsive to 
climate change. Increased warming during the growing season 
is likely to alter production by lengthening the growing season 
and increasing metabolic activity [3]-[5]. Changes in seasonal 
air temperature may increase soil temperatures and the active 
layer depth to permafrost. This, in turn, affects soil microbial 
activity, nutrient cycles and soil moisture [6], [7], which alter 
existing plant growth and competitive species interactions re- 
sulting in dramatic changes in vegetation composition. Further, 
tundra ecosystems are often moisture limited, so climate in- 
duced changes in precipitation patterns and surface hydrology 
will also act to alter vegetation growth patterns [8], [9]. 

A unique characteristic of the tundra ecosystem is the rela- 
tive dominance of non-vascular plants in the landscape, such 
as mosses and lichens. In some cases, mosses dominate local 
habitats in the tundra because of their tolerance to extreme cold, 
dessication, water-logging, and low light [10]. Lichens can sur- 
vive extreme climate conditions such as drought and persistent 
freezing as well [1 1]— [13]. Tundra warming manipulation ex- 
periments have shown a decline in lichens with an increase in 
vascular plants suggesting changing proportions of tundra plant 
types as climate changes [14]. 

Non-vascular plants can provide a significant fraction of 
tundra carbon uptake and should be explicitly included in 
descriptions of tundra carbon fluxes. Photosynthesis of tundra 
vascular plant canopies has been described with strong rela- 
tionships found between leaf area or the fraction of absorbed 
photosynthetically active radiation and gross primary produc- 
tivity [15], [16], [3]. However, the net primary productivity 
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(NPP) of mosses have been shown to represent about 25-30% 
of total above-ground NPP in several tundra sites including; 
tussock tundra [9], coastal tundra [17] and tundra heath [18]. In 
other cases mosses were found to dominate tundra photosyn- 
thetic C0 2 fluxes [19], [20]. 

There are further important physiological differences be- 
tween vascular and non- vascular plants that affect tundra water 
and energy balance. Mosses and lichens do not have roots, so 
cannot access soil moisture below the surface. Instead they 
depend on water from atmospheric humidity, ground water, or 
precipitation. They also do not possess stomata and therefore 
lose water readily to the atmosphere and many species are 
adapted to survive long periods of desiccation. [21]. They may 
also act as an insulating layer to heat transport between the soil 
and the atmosphere [22]. Thus, in areas where these non-vas- 
cular plants make up a significant portion of the vegetation 
cover the response of ecosystem carbon and energy balance 
to environmental changes is expected to differ from that of 
vascular plant-dominated vegetation and must be accounted for 
in tundra ecosystem modeling. 

The difficulties in working in remote tundra locations have 
always incurred large logistical costs, therefore observations 
of ecosystem change using remote sensing provide an expe- 
dient and economical method of collecting repeatable and 
consistent measurements over large areas. Importantly, most 
remote sensing theory of vegetation has been derived from 
studies focused on temperate crops and forests and thus the 
interpretation of conventional remote sensing tools is often 
ambiguous for northern landscapes. There are several unique 
features of the tundra that affect the interpretation of remotely 
sensed data [23]. Non-vascular plants can represent significant 
fractions of the tundra landscape cover, and these non-vascular 
plants have different spectral characteristics as well as different 
physiological responses from vascular plants or from bare soils 
[24]-[26]. For example, remote sensing studies have been able 
to detect changes in lichen-dominated areas due to responses 
to both short-term temperature anomalies as well as long-term 
temperature trends [27], [28]. 

This study investigates the importance of non-vascular plants 
in the tundra landscape and examines their effects on tundra 
carbon uptake, using spectral reflectance to distinguish these 
different functional types. We explore the concept of “optical 
types” [29] to scale from ground-based in situ measurements to 
landscapes using satellite observations to examine vegetation 
functional type patterns across the tundra. 

II. Materialsand Methods 

A. Study Area 

The study examines an area mainly located north and west of 
the town of Barrow, AK, including most of the Barrow Envi- 
ronmental Observatory. Specific study plots are located at 71° 
19’ 19.1” N 156° 36’ 15.9” W, approximately seven kilome- 
ters east of Barrow. The area is classified as moist acidic coastal 
tundra [30]. Barrow has a mean annual temperature of — 12.0°C, 
with a minimum mean monthly temperature of — 26.6° C in Feb- 
ruary and maximum mean monthly temperature of 4.7° C in 
July. More than half of the 1 06 mm of annual precipitation falls 


as rain during a 3 -month period from July through September, 
with the ground being snow covered from September through 
May to mid- June [31]. 

The tundra community at Barrow, AK consists of an over- 
story dominated by vascular plants (mostly graminoids) and 
an underlying mat of mosses [32]. Vascular plants at Barrow 
include dwarf shrubs, forbs, and perennial herbs but the most 
dominant growth form is graminoid [30]. Common vascular 
species are: Carex aquatilis, Dupontia fischerii, Eriophorum 
angustifolium, E. scheuzerii, Luzula confusa, Petasites frigidus, 
Potentilla hyparctica, and Salix rotundifolia. Mosses are a sig- 
nificant portion of the community and in some areas may 
account for more biomass than vascular plants [33]. Moss 
species include: Dicranum elongatum, D. undulatum, Drepan- 
ocladus revolvens, Polytrichum (juniperum), and Sphagnum 
spp. Lichens are scattered throughout the drier areas but are not 
as abundant or productive as mosses and graminoids [34] and 
include the species: Alectoria nigricans, Cetraria cucullata, C. 
nivalis , and Dactylina arctica. 

B. Measurements 

This project included: 1) the collection of in situ mea- 
surements of plots for the three vegetation functional types, 
including photosynthesis and spectral reflectance; 2) measure- 
ments of ground cover amount and spectral reflectance at every 
meter along a 100 m transect; 3) development of relationships 
between spectral reflectance and both photosynthesis and cover 
amount; 4) estimation of LUE across a tundra transect; and 
5) application of these findings to the spectral information 
collected from Hyperion. 

Vegetation measurements were designed to relate photo- 
synthetic rate to spectral reflectance. Tundra vegetation were 
grouped into three functional groups; vascular plants, mosses, 
and lichens, with at least four replicates per group measured on 
any one day. Each replicate provided a different species from 
each group in an attempt to measure the variability within a 
group. Field sampling for this study was conducted on July 20 
and August 5, 2001, near the peak of the growing season in this 
area. 

Spectral reflectance measurements were collected at 5 1 plots, 
each about 15 cm in diameter, representing nearly pure areas 
of representative tundra species. Photo synthetic gas exchange 
measurements were collected at 17 of the plots. All of the study 
plots were in very close proximity (< 15 m) to a 100 m linear 
transect that was also optically sampled at frequent intervals 
throughout the summer [35]. 

CO 2 gas exchange was measured at 17 plots where a 14.6 cm 
diameter collar of acrylic plastic was embedded in the ground. 
When measurements were made, the plot was enclosed in a clear 
acrylic chamber constructed of 0.3175 cm thick acrylic tube 
glued to a 0.635 cm thick acrylic top and a closed gas exchange 
system was connected to the chamber (LI-6200, LiCOR, Lin- 
coln, NE). Each individual plot was measured under ambient en- 
vironmental conditions for net photosynthesis and then the mea- 
surement was repeated with the chamber covered with a black 
cloth to block out light to provide an estimate of dark respira- 
tion. Gross photosynthesis was calculated as the difference be- 
tween the net photosynthesis, measured with the clear chamber, 
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and dark respiration, measured with the chamber covered. A 
quantum sensor mounted on the chamber top measured incident 
photosynthetically active radiation (PAR). The chamber had a 
PAR transmittance of 95% [36]. Photo synthetic light-response 
curves were developed using household window screen material 
as neutral density filters to control incident PAR, allowing pho- 
tosynthesis measurements for each functional type to be rapidly 
collected under a range of light levels. 

Light use efficiency (LUE) was determined as the ratio of 
absorbed light and gross photosynthesis: 


/par Qi 

where e is LUE, G is gross photosynthesis with units of 
//mol CO 2 m -2 s -1 , Qi is incident PAR with units of //mol 
photons m -2 s -1 (i.e., photo synthetic photon flux density, or 
PPFD) and /par is the fraction of the PAR absorbed by the 
vegetation (unitless). Absorbed PAR (APAR) is the product of 
/par and incident PAR. 

/par is not well defined and difficult to determine for vege- 
tation with low growth forms such as lichens and mosses. We 
assumed that no PAR was transmitted through moss and lichen 
mats, so for these vegetation types /par is 1-PAR albedo. PAR 
albedo is estimated by integrating the measured spectral re- 
flectance (see below) between 400 and 700 nm. For the vas- 
cular plant plots /par is determined using the linear relation- 
ship between Normalized Difference Vegetation Index (NDVI) 
and /par described in Huemmrich et al [36]. 

The optical field sampling consisted of measurements of the 
small single species plots, as well as measurements collected at 
every meter along a 100 m transect to observe areas of mixed 
coverage [35], [37], [38]. Optical measurements were collected 
using two portable field spectrometers (UniSpec, PP Systems, 
Haverhill, MA). These spectrometers measured reflectance be- 
tween 310 and 1130 nm sampling at approximately 3 nm inter- 
vals. 

For plot measurements, the spectrometer was fitted with a 
glass fiber optic cable connected to a stainless steel ferrule, 
which provided a 20° field of view. Reflected irradiance mea- 
surements were collected with the end of the fiber optic cable 
held vertically between 30 and 40 cm above the surface to view 
an area of less than 1 5 cm diameter to match the area of gas ex- 
change measurements. Within seconds of collecting a surface 
measurement, a reference measurement was made of a calibra- 
tion panel (Spectralon, Labsphere, Inc. North Sutton, NH). The 
speed of this procedure allowed the acquisition of data even 
under cloudy conditions. Three reflectance spectra were col- 
lected for each sampled plot, and were averaged to produce 
one spectrum per plot. In addition to the vegetation plots, spec- 
tral measurements for nongreen materials including bare soil, 
open water, and areas of standing dead vascular plants, were 
collected soon after snowmelt and before the start of green-up. 
For the vegetation plots, ground-based spectral reflectance mea- 
surements were collected near the height of peak productivity 
period in early August (August 5 and 8, 2001) and reflectance 
for nongreen plots were collected throughout the summer (June 
20 and 30, July 15, August 8 and 22, and September 3, 2001). 


To characterize spectral reflectance patterns at multiple spa- 
tial scales, spectral measurements were made at 1 m intervals 
along a 100 m linear track system [35], [38]. The track was sup- 
ported less than a meter above the tundra surface by tripods with 
a tram cart riding on the track carrying a dual channel spectro- 
radiometer (Unispec DC, PP Systems). From the spectrometer 
were two fiber-optic cables mounted on a mast and boom at- 
tached to the cart, one attached to a diffuser head viewing up- 
ward, the other on the boom pointing vertically downward to 
view an area south of the track with a field-of-view of approxi- 
mately 1 m in diameter. With each measurement of reflected ra- 
diance, a measurement of incident irradiance was also acquired 
at the same time for a determination of surface reflectance. Due 
to the coincident measurements of both incoming and reflected 
radiation this system was able to correct for effects of varying 
cloud cover. The tram system collected spectral measurements 
that were repeatable throughout the growing season while mini- 
mizing disturbance to the tundra surface [35], [38]. In this study 
we used measurements from the peak of the growing season pe- 
riod collected on August 10, 2001. 

Processing the ground-measured spectral reflectance data in- 
cluded interpolating the spectral bands to 1 nm intervals. Due 
to instrument noise at the ends of the spectrometer range, the 
usable spectral range was between 400 and 1 000 nm. The pro- 
cessing used freely available software (Multispec v.4.0, http:// 
specnet.info/specnet_toolkit.htm). The spectral reflectance data 
were then convolved to match the approximately 1 1 nm band 
passes of the Hyperion instrument [39] in subsequent analyses. 

Imaging spectrometer data were acquired on July 20, 2009 
(day 20 1 ) by the Hyperion instrument on the Earth Observing- 1 . 
Hyperion provides data in 220 spectral bands from 0.4 to 2.5 //m 
at 30-meter resolution and a 7.5 km wide swath [39]. The satel- 
lite data were atmospherically corrected to surface reflectance 
using the Atmosphere Removal algorithm (ATREM) [40], [41]. 
To match the spectral range of the ground-based spectral data 
only 55 of the continuous Hyperion bands between 437 and 993 
nm were used. 

Discriminant analysis was used to separate the functional 
groups based on optical properties, determining the important 
spectral bands and simplifying the number of independent 
variables [42]. The plot data were divided into four groups: 
vascular plants, mosses, lichens, and nongreen materials (bare 
soil, dead vegetation, and standing water). Inputs to the dis- 
criminant analysis using Systat (Version 13, Systat Software, 
Inc., Chicago, IL) were plot reflectance spectra convolved to 
the Hyperion bands. The analysis used a forward stepwise 
method (Wilks’ Lamda) to test for significant mean group 
differences, using a probability of F < 0.05 to include bands 
and F > 0.10 to remove the band. 

At each meter along the 100 m transect a visual estimate of 
species abundance was made over the period August 8-11, 200 1 
using a 1 m by 1 m quadrat. The species data were converted 
into percent area of each of the vegetation functional types and 
nongreen material types. We estimate errors in the estimation of 
coverage to be ±15% cover. Because of the three dimensional 
structure of the vegetation the sum of the coverage was more 
than 100%. For this analysis, the total measured cover values 
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were normalized so total cover equals 1 00% as is generally done 
in linear unmixing approaches. 

LUE was determined for areas with mixtures of different 
functional types as the weighted sum of the LUE from the pure 
plots, where the weight coefficients were the normalized frac- 
tional cover values. This calculation assumes that the incident 
PAR is absorbed by the landscape components proportional to 
their normalized cover fraction and was applied to determine 
LUE for both the transect data and the Hyperion satellite 
imagery. 

Chlorophyll concentrations are a key physiological factor re- 
lated to maximum photo synthetic rates. We estimated chloro- 
phyll concentrations using Gitelson’s three-band model: 

Ci=(- M PA3 (2) 

\p\l p\2j 

where the spectral index Q is proportional to the chlorophyll 
concentration, and p\ is the reflectance for a given wavelength 
band, A. Wavelengths chosen are maximally sensitive to ab- 
sorption by chlorophyll and other pigments (A1 = 549 nm), 
are influenced primarily by non-chlorophyll pigments (A2 = 
793 nm), and where reflectance is controlled by leaf scattering 
(A3 = 793 nm) [43]-[45]. 

The cover estimates of the functional types retrieved from the 
Hyperion imagery were used to derive landscape LUE patterns. 
The LUE per pixel was calculated as described above and com- 
pared with chlorophyll concentration estimates obtained using 
( 2 ). 

III. Results 

The physiological measurements identify differences re- 
lating to photosynthetic carbon exchange among the three 
vegetation functional types. From the plot gas exchange data, 
LUE for each functional type was calculated as the slope of 
a linear regression between gross photosynthesis and APAR, 
forced through the origin (Fig. 1, Table I). Lichens, mosses, and 
vascular plants all had different LUE values. Vascular plants 
had significantly higher net photo synthetic rates and LUE 
than the other two functional types. The vascular plant LUE 
from this study (0.0134 mol C mol -1 quanta) was close to 
that previously determined for a vascular plant overstory near 
Barrow (LUE = 0.0126 mol C mol -1 quanta [36]) and within 
the 95% confidence interval of the slope from the regression in 
this study (i.e., the vascular plant LUE). This was noteworthy 
as the LUE in this study was based only on measurements 
collected during the middle of the growing season, while the 
Huemmrich et al. [36] LUE was derived from measurements 
collected throughout the growing season, suggesting the LUE 
for vascular plants was relatively stable over much of the 
season, as well as having similar LUE values for vascular 
plants in mixtures with mosses or growing alone. 

The LUE for mosses and lichens were much lower than for 
vascular plants. There was a significant difference, however, be- 
tween these two functional types (Table I) with mosses having 
a LUE almost 50% greater than lichens (a similar analysis for 
LUE of the moss data was presented in [36]). 

The in situ spectral reflectance measurements of the plots 
were divided up into four groups: vascular plants (n = 22), 



Absorbed PAR (pmol nr 2 s' 1 ) 


Fig. 1 . Absorbed PAR and gross photosynthesis for the three vegetation func- 
tional types from plot data collected July 20 and August 5, 2001. Light use 
efficiency for each functional type was calculated from these data as the slope 
forced through the origin. See Table I for description of regressions. 


TABLE I 

Light Use Efficiency Based on Regressions Between Net 
Photosynthesis and APAR From Plot Data. LUE is in Units of 
mol C mol -1 Absorbed Quanta 



Vascular 

Plants 

Moss 

Lichen 

LUE 

0.01337 

0 00184 

0.00124 

St. Err, of 
Slope 

0.00101 

0,00022 

0.0001 1 

Num. Obs. 

29 

25 

28 

R 3 

0.74 

0.56 

0.74 

P 

<0.01 

<0.01 

<0.01 


mosses (n = 11), lichens (n = 18), and nongreen material 
(n = 19) (Fig. 2). Generally, the vascular plant and moss 
spectra were similar, and display typical green plant spectral re- 
flectance patterns, including more variation in the near infrared 
than in visible wavelengths. In comparison, lichen spectra had 
much higher visible reflectance and were more variable in 
all wavelengths. Spectral patterns for the nongreen materials 
were clearly different for dead vegetation, bare soil, and water. 
Generally the spectral variability for the entire nongreen group 
increases as a function of wavelength. 

Discriminate analysis provided an objective tool to evaluate 
the ability of optical sampling to distinguish functional groups 
and to identify a subset of the spectral bands to do this separa- 
tion. The stepwise method for discriminate analysis reduced the 
number of wavebands from 55 to five, yet retained separation 
of functional types. The five chosen bands were located in key 
locations of typical green plant spectra (Fig. 2); at the blue (488 
nm) and red (671 nm) chlorophyll absorption wells, near the in- 
flection point of the red edge (712 nm), at the shoulder of the 
red edge (763 nm), and in the near infrared “plateau” (834 nm). 
Mulhern [46] identified the blue band as an important spectral 
region for separating lichens from soils. The discriminate anal- 
ysis functions (Fig. 3) correctly predicted 83% of group mem- 
bership overall, with 86% correct for vascular plants, 91% for 


HUEMMRICH et al . : ARCTIC TUNDRA VEGETATION FUNCTIONAL TYPES BASED ON PHOTOSYNTHETIC PHYSIOLOGY AND OPTICAL PROPERTIES 


269 



Wavelength (nm) 

Fig. 2. Endmember spectral reflectance for each cover type, vertical lines indicate wavelengths used in the discriminant analysis. 



Discriminant Function 1 

Fig. 3. The distribution of data from plots using the first two functions derived 
from discriminate analysis for the three vegetation functional types and non- 
green materials. The larger gray circles indicate values used as endmembers in 
the unmixing analysis. For the endmember points the symbol within the circle 
indicates the vegetation functional type of that endmember. 

mosses, 95% for the nongreen materials, and 61% for lichens. 
Lichen errors were mainly due to confusion with the nongreen 
materials (Fig. 3). 

Even at the scale of one square meter, the cover estimates 
along the tram transect show that all square meter plots 
consisted of mixtures of multiple cover types (Fig. 4). The 
observed cover estimates included little bare ground with max- 
imum coverage of 20%, while standing dead vegetation ranged 
from 0-70%, and water from 0-100%. Moss coverage ranged 
from 0-100%, but when normalized to a 100% maximum total 
coverage maximum moss coverage was 40% 4= 7%, based on 
the propagation of an assumed 15% error in the original obser- 
vations. Lichen coverage ranged from 0-90% with 35% ± 7% 
maximum normalized lichen coverage. Vascular plant cov- 
erage ranged between 40-150%, with normalized coverage 
of 18-88% zb 9%. Vascular plant coverage was highest in the 
locally low areas, particularly the areas with standing water, 
while moss coverage was highest in locally low areas without 
standing water and on the edges of the wet areas. Lichen 
coverage was highest in the locally high areas where bare soil 


patches also occurred. These spatial patterns indicate the role 
of microtopography on the distribution of the functional types. 

The transect coverage observations indicate how heteroge- 
neous the tundra is, as areas over a few square centimeters in size 
generally consist of mixtures of multiple types. Therefore, to es- 
timate LUE at that scale of a 30 m Hyperion pixel, the fractions 
of each functional type must be determined. We used the op- 
tical measurements to scale from plot-level LUE to 30 m areas. 
Since the discriminant functions optimized separation of the dif- 
ferent functional types with a small number of spectral bands, 
we used the distance to endmembers in the discriminate anal- 
ysis function space to estimate cover fractions in mixed pixels. 
Assuming that the range of in situ plot reflectances provide a 
reasonable description of the variability of plot type, the end- 
members for the moss, lichen, and nongreen cover types were 
assigned values of the average of the plot reflectances. However, 
variability in the vascular plants reflectance is driven by differ- 
ences in green leaf area index. To account for that variability the 
vascular plant endmember was chosen to be the “greenest” plot 
spectra (i.e., spectra with the highest NDVI, which was also the 
greatest distance from the average of all the plot values in the 
discriminate analysis function space). 

Endmember reflectance for each functional type and re- 
flectances at each meter along the transect were transformed 
using the discriminate analysis functions. The distances in the 
discriminate analysis space between the coverage estimates 
for each square meter and the endmembers were calculated, 
these were related to the observed coverage to create transfer 
functions. The statistical associations between distance from 
endmembers and observed coverage were low due to large 
observation errors relative to the range of values, with R 2 
values of 0.39, 0.13, and 0.19 for vascular plants, mosses, 
and lichens, respectively (Table II). Nevertheless, the cover 
fractions determined from the in situ reflectance were found to 
be in close agreement with sampled cover fractions (Fig. 6), 
with root mean square error (RMSE) of 11% cover for vascular 
plants, 8% for moss, 9% for lichen, and an overall RMSE of 
9% cover. 

LUE was calculated for each square meter block along the 
transect using the observed cover fraction and coverage esti- 
mated from the reflectance spectra (Figs. 5 and 7). Along the 


270 


IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING, VOL. 6, NO. 2, APRIL 2013 



■ Tola I VP 
■Moss 
Uchen 
Nongneen 
■Water 


Fig. 4. Observed functional type coverage at 1 m intervals along the transect, visually estimated over the period August 8-11, 2001. Because the vegetation is 
3-dimentional, total cover fractions were over 100%. In the data analysis the total cover fractions were normalized to equal 100%. 
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Fig. 5. LUE along the 100 m transect shown with the black line was estimated from the functional type coverage (shown in Fig. 4) and light use efficiency (shown 
in Fig. 1), along with LUE estimated using coverage derived from spectral unmixing shown with the gray line. Microtopography is shown as dashed line. Error in 
Observed LUE estimate is 0.00136 mol C mol -1 quanta, RMSE in LUE derived from observations compared to LUE from spectra is 0.00266 mol C mol -1 
quanta. 


TABLE II 

Description of Relationships in Spectral Space Between the 
Statistical Distance to Endmembers and Coverage Estimate 
Observations, Acquired Along the Tram Transect. The Errors 
are Expressed in Percent Coverage. For all Cases the Number of 
Observations is 100 



Vascular 

Plants 

Moss 

Lichen 

R* 

0,39 

0.13 

0,19 

P 

<0.01 

<0.01 

<0.01 

Esl. Obs. Error 

9,5 

7.1 

7.0 

SE of Regression 
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transect LUE was quite variable, with maximum and minimum 
values occurring within 6 m of each other (Fig. 7). Variability 
in LUE was affected by microtopography, with higher values of 


LUE tending to occur in locally low, wet areas. This variability 
over the 1 00 m distance produced a fourfold difference in LUE 
with values from 0.003 to 0.012 mol C mol -1 quanta for the in 
situ observations. The average errors for estimating LUE were 
similar for the two methods: 0.0014 mol C mol -1 quanta for 
the method using observed coverage and 0.0015 mol C mol -1 
quanta for the LUE based on the remote sensing approach. The 
two different LUE calculations have a correlation of 0.62 with 
a RMSE of 0.0013 mol C mol -1 quanta. 

The equations relating spectral reflectance and functional 
type coverage derived from the transect data were the applied to 
the Hyperion imagery to create a continuous fields description 
of the coverage of the functional types for the tundra around 
Barrow (Fig. 8(b)). Over this area the distribution of vascular 
plant coverage per pixel had a maximum value of 74% with a 
distribution peak of 32%, for mosses the maximum coverage 
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Cover Fraction from Spectra 

Fig. 6. Comparison between observed vegetation functional type coverage nor- 
malized to 100% for vascular plants, mosses, and lichens along the transect and 
cover fractions estimated from spectral reflectance. Line is 1 to 1 line. The ab- 
solute error in visual estimates is approximately 9%. The RMSE of spectral 
retrievals is 9% absolute. 



LUE from Spectral Cover Estimates 

Fig. 7. LUE (mol C mol - 1 quanta) along the transect estimated based on ob- 
served normalized cover fractions and cover fractions obtained from spectral 
unmixing. Line is 1 to 1 line. Errors in LUE based on observed cover fractions 
are due to errors in LUE for the pure functional types and an assumed 15% error 
of observed raw coverage amounts. Errors in LUE based on spectrally derived 
cover fractions are due to errors in LUE for pure functional types and the RMSE 
of the spectral retrieval of cover fractions. 


was 31% with the distribution peak at 22%, and for lichens the 
maximum coverage was 22% with a peak of 10%. 

Without ground observations of vegetation type coverage, an 
examination of an area known to be free of vegetation was used 
to make an error estimate. The old Naval Arctic Research Labo- 
ratory runway was chosen as for this test. For runway pixels vas- 
cular plant coverage was as high as 9%, moss 14%, and lichen 
13%, compared to expected values of zero percent coverage. 


These values were close to the RMSE of coverage from the tram 
analysis comparing the spectral unmixing with ground obser- 
vations, suggesting the Hyperion unmixing captured the actual 
functional type distribution near this level of error. 

There are clearly observed spatial patterns in the cover frac- 
tions in this landscape. High vascular plant coverage occurred 
on the margins of drained lakes. Higher lichen coverage tended 
to be found in drier upland regions. This pattern of lichen 
occurrence was also observed in the transect measurements 
(Fig. 4). Moss coverage was more widely distributed, although 
high moss coverage in the middle of the drained lakes may be 
erroneous due to confusion between mosses and mixtures of 
water and vascular plants. Field observations confirm that moss 
often occurs as a low, “background” layer beneath the vascular 
canopy of wet tundra in this region [36]. 

The Hyperion functional type coverage estimates were used 
to calculate landscape patterns of LUE (Fig. 8(c)). For the veg- 
etated parts of the scene, this Hyperion-based LUE ranged from 
0.0021-0.0102 mol C mol -1 quanta, a fivefold difference be- 
tween maximum and minimum LUE. The distribution peak of 
LUE for this area and the overall average were both 0.0048 
mol C mol -1 quanta. 

The Hyperion-based LUE was compared with the chlorophyll 
spectral index, Ci ((2)). Ci was well correlated with LUE over 
the study area (R 2 = 0.69, SE = 0.0006 mol C mol -1 quanta) 
(Fig. 9) suggesting vegetation chlorophyll concentration is a 
key determinant of LUE for this tundra ecosystem. This sug- 
gests that a more direct approach for deriving spatial patterns of 
LUE for tundra from imaging spectrometry would be based on 
a simple reflectance-based metric of chlorophyll concentration. 

IV. Discussion and Conclusions 

This study uses the concept of optical types [29] to examine a 
key characteristic of ecosystem carbon exchange in the tundra. 
To have usable optical types one must be able to both identify 
significant functional differences in vegetation types and spec- 
trally identify these different types. Our grouping of tundra veg- 
etation (vascular plants, mosses, and lichens) meets these cri- 
teria. Among the different types, the plot measurements showed 
distinct differences in LUE, a key variable describing photosyn- 
thetic carbon uptake and discriminate analysis showed that the 
types could be separated based on their spectral reflectance. The 
ability to link LUE to reflectance characteristics provides an ap- 
proach to scale from the ground measurements to distributions 
over a landscape. Using an unmixing approach we were able to 
create continuous fields of key tundra functional types for this 
area, an improvement over simple classifications. Results from 
this study demonstrated how, even at a scale of a square meter, 
variable mixtures of functional types produced significantly dif- 
ferent values for LUE. Unmixing using spectral reflectance pro- 
vides a unique estimate of LUE for each pixel. Remote sensing 
provides a way to map large areas and make noninvasive repeat 
measurements to monitor ecosystem change. Remote sensing 
is a particularly important tool for observing tundra due to the 
difficulties and expense involved in traveling to and working in 
this region, rendering direct field sampling over large areas im- 
practical. 
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Fig. 8. Three versions of the EO-1 Hyperion image acquired on July 20, 2009 are shown. Left: 3-band (RGB = 834, 671, and 549 nm) composite image of 
surface reflectance. The grid of light blue lines on the lower left side of the image is the city of Barrow. The straight blue line along the shore near the top of the 
image is the old airport runway, used in the error evaluation. The oblong features scattered around the region are drained thermokarst lakes and the dark red ones are 
now marshes. Middle: Three band RGB continuous fields of estimated coverage of vegetation functional types derived from spectral unmixing and scaled between 
0 and 50% coverage. R = Vascular Plant Cover, G = Moss Cover, B = Lichen Cover. Right image: Map of LUE spatial patterns (mol C mol -1 quanta 
xlOOO) based on coverage estimates. 
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Fig. 9. Scatter plot of chlorophyll index (x-axis) versus LUE (mol C mol -1 
quanta xlOOO) based on normalized cover fractions (y-axis). Both variables 
were retrieved from the Hyperion image. 


The linear unmixing approach used to calculate pixel-level 
LUE simplifies the tundra canopy as it does not take into ac- 
count the three dimensional structure of a vascular plant over- 
story. Vascular plants can grow through mats of mosses and 
lichens and shade them, decreasing their photosynthetic pro- 
duction [36]. However, vascular plant leaf area index for this 
region is generally low, with observed midseason values less 
than 2 [36]. We believe the overall effect of this overstory on 
the moss and lichen understory production is relatively small in 


this region due to both the low leaf area combined with gener- 
ally diffuse light conditions. 

This study provided a snap-shot of this landscape at a par- 
ticular point in time. Seasonal growth of vascular plants will 
change their relative proportions throughout the growing season 
and this affects carbon uptake and LUE patterns [38]. Early in 
the growing season, as vascular plants begin to grow, tundra pro- 
ductivity will be dominated by moss and lichen photosynthesis, 
with the vascular plant component increasing as the season pro- 
gresses [36]. Thus, the spatial patterns of coverage and LUE are 
expected to change through the growing season. 

The correlation between LUE calculated from coverage 
estimates and the chlorophyll index (Ci), even though they 
use different algorithms and spectral bands, demonstrates how 
multiple approaches can yield similar results. This convergence 
is most likely due to plant physiological constraints indicating 
chlorophyll concentration is a controlling factor in vegetation 
photosynthetic efficiency, and may provide a more direct 
approach to deriving LUE. 

This study found significant variability in functional type 
cover fractions, leading to variability in estimated LUE, at local 
scales (meters) in the field data collected along the transect and 
at landscape scales (km) in the satellite imagery. The spatial 
variability shown in LUE estimates is not accounted for in 
existing carbon flux models. Spatial distributions of vegetation 
functional types were strongly related to surface microtopog- 
raphy, from which we infer surface hydrology to be a significant 
controlling factor. These results suggest that climate change 
affecting surface hydrology [47] will also affect spatial patterns 
of vegetation distribution and ecosystem carbon exchange. 
Particularly in a scenario of modified surface hydrology, we 
can expect that climate change will affect relative coverage due 
to differing responses of each of the functional types to climate 
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change, and that these changes in cover will be associated with 
altered carbon and energy balance. 

Since hydrology is an important factor determining the vege- 
tation functional type coverage, shortwave infrared bands may 
provide additional information on vegetation water content, im- 
proving the retrieval of functional type coverage [48]. This ap- 
proach was not tested in this study as the field measurements 
observed only visible and near infrared spectral bands. Recent 
observations combining NIR and blue bands to detect standing 
water in coastal tundra [49] may also prove useful at detecting 
hydrological changes linked to cover type and LUE. 

The results of this study are indicative only of this part of 
the Alaskan coastal tundra and further field studies are required 
in other parts of the tundra biome to test and extend the re- 
sults and to evaluate the ability to generalize the results. For ex- 
ample, shrubs are an important component of the tundra biome 
and shrub coverage has been shown to be increasing in the 
tundra [50]. Shrubs did not represent a significant cover type 
in our study area, however, so were not explicitly included in 
this study. Future work is required to examine shrub optical and 
physiological characteristics to see if they comprise an optically 
distinct functional type. 

The use of multiple narrow spectral bands for the umixing 
points to an application for the future NASA Hyperspectral In- 
frared Imager (HyspIRI) mission for mapping tundra vegetation 
distribution according to physiological function, hydrology, and 
micro topography over the entire biome. The concept of optical 
types (optically distinguishable functional types) may prove to 
be useful for this kind of functional mapping. Similarly, simple 
metrics of surface water cover (e.g., [49]) or pigment distribu- 
tion (e.g., [45]) may provide insight into the proximal causes of 
functional change in tundra ecosystems. 

The close linking of optical properties and carbon flux mea- 
surements in this study has provided insights into the func- 
tioning of the tundra ecosystem by identifying fine-scale pat- 
terns of FUE that may be related to hydrology and microto- 
pography. This study also illustrates new “scaleable” functional 
mapping techniques that can be extended through the use of 
satellite remote sensing to larger regions. This approach lends 
itself to future studies considering how species composition re- 
lates to changing ecosystem function and provides a framework 
for studying ecosystem change through remote sensing that con- 
siders shifts in hydrology, species composition, and their effects 
on carbon balance. 
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